function [zeta,weight]=GetTriGaussIntgInfo(GsNum,inp)
zeta=zeros(2,1);
if GsNum==3
    coords=[0.66666666666666666666666666666666666667, 0.16666666666666666666666666666666666667;
        0.16666666666666666666666666666666666667, 0.66666666666666666666666666666666666667;
        0.16666666666666666666666666666666666667, 0.16666666666666666666666666666666666667];
	weights=[0.3333333333333333, 0.3333333333333333, 0.3333333333333333];
    zeta(1)=coords(inp,1);
    zeta(2)=coords(inp,2);
    weight=weights(inp);
elseif GsNum==4
    coords=[0.6, 0.2;
        0.2, 0.6;
        0.2, 0.2;
        0.333333333333333, 0.333333333333333];
	weights=[0.5208333333333333, 0.5208333333333333, 0.5208333333333333, -0.5625];
    zeta(1)=coords(inp,1);
    zeta(2)=coords(inp,2);
    weight=weights(inp);
elseif GsNum==7
    coords=[0.333333333333333, 0.333333333333333;
			0.101286507323456, 0.101286507323456;
			0.797426985353087, 0.101286507323456;
			0.101286507323456, 0.797426985353087;
			0.470142064105115, 0.470142064105115;
			0.059715871789770, 0.470142064105115;
			0.470142064105115, 0.059715871789770];
        weights=[0.225000000000000
				0.125939180544827
				0.125939180544827
				0.125939180544827
				0.132394152788506
				0.132394152788506
				0.132394152788506];
    zeta(1)=coords(inp,1);
    zeta(2)=coords(inp,2);
    weight=weights(inp);
else
    
end
        
